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Samples of inclusive 7+2 jet and 7+3 jet events collected by the DO experiment with an integrated 
luminosity of about 1 fb~^ in pp collisions at = 1.96 TeV are used to measure cross sections as a 
function of the angle in the plane transverse to the beam direction between the transverse momentum 
(Pt) of the 7+leading jet system (jets are ordered in pr) and pr of the other jet for 7 + 2 jet, or 
Pt sum of the two other jets for 7 + 3 jet events. The results are compared to different models of 
multiple parton interactions (MPI) in the pythia and SHERPA Monte Carlo (MC) generators. The 
data indicate a contribution from events with double parton (DP) interactions and are well described 
by predictions provided by the pythia MPI models with pT-ordered showers and by SHERPA with 
the default MPI model. The 7 + 2 jet data are also used to determine the fraction of events with 
DP interactions as a function of the azimuthal angle and as a function of the second jet pr- 

PACS numbers: 13.85.Qk, 12.38.Qk 



I. INTRODUCTION 



The high energy scattering of two nucleons can be con- 
sidered, in a simplified model, as a single collision of one 
parton (quark or gluon) from one nucleon with one par- 
ton from the other nucleon. In this approach, the remain- 
ing "spectator" partons, which do not take part in the 
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hard 2 — 2 parton collision, participate in the so-called 
"underlying event." However, there are also models that 
allow for the possibility that two or more parton pairs un- 
dergo a hard interaction when two hadrons collide. These 
MPI events have been examined in many theoretical pa- 
pers (l|-[T7j. A comprehensive review of MPI models in 
hadron collisions is given in Q. A significant amount of 
experimental data, from the CERN ISR pp collider [l^ , 
the CERN SPS pp colhder [11], the Fermilab Tevatronnp 
coUider [13^ and the DESY HERA ep colhder [H,!!!, 
shows clear evidence for MPI events. 

In addition to parton distribution functions (PDF) and 
parton cross sections, the rates of double and triple par- 
ton scattering also depend on how the partons are spa- 
tially distributed within the hadron. The spatial parton 
distributions are implemented in various phenomcnologi- 
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cal models that have been proposed over the last 25 years. 
They have evolved from the first "simple" model sug- 
gested in [3] , to more sophisticated models [HI that 
consider MPI with correlations in parton momentum and 
color, as well as effects balancing MPI and initial and 
final state radiation (ISR and FSR) effects, which are 
implemented in the recent ( "pT-ordcred" ) models [T2j . 

Beyond the motivation of better understanding non- 
perturbative quantum chromodynamics (QCD), a more 
realistic modeling of the underlying event and an esti- 
mate of the contributions from DP interactions are im- 
portant for studying background events for many rare 
processes, including searches for the Higgs boson 113- 
l3l| . Uncertainties in the choice of the underlying event 
model and related corrections also cause uncertainty in 
measurements of the top quark mass. This uncertainty 
can be as large as 1.0 GeV [s^, a value obtained from 
a comparison of MPI models with virtuality-ordered @ 
("old" models) and pT-ordercd [l^ ("new" models) par- 
ton showers. 

In a previous paper p^ . we have studied the 7-1-3 jet 
final state and extracted the fractions of DP events from 
a comparison of angular distributions in data with tem- 
plates obtained from a data-driven DP model. That pa- 
per also presented the effective cross section ((Toft), which 
characterizes the size of the effective parton-parton inter- 
action region and can be used to calculate the DP cross 
sections for various pairs of parton scattering processes. 

In this paper, we extend our previous study by mea- 
suring differential cross sections for the angle in the plane 
transverse to the beam direction between the pt vector 
obtained by pairing the photon and the leading (ordered 
in Pt) jet and the pr vector of the other (two) jct(s) 
in 7 + 2(3) jet + X events (referred to below as "7 + 
2(3) jet" events). These cross sections are sensitive to 
the contributions from jets originating from additional 
parton hard interactions (beyond the dominant one) and 
can be used to improve existing MPI models, and to es- 
timate the fractions of such events [2l|, [13, . The cross 
section measurements are performed at the particle level, 
which means that the jets' four-momenta represent the 
real energy and direction of the jet of stable particles re- 
sulting from the hadronization process following the pp 
interaction [s^l ■ The larger statistics in 7 -I- 2 jet events 
allows to subdivide the cross section measurement in bins 
of the second jet pT (p^^*'^). This extension increases the 
sensitivity to various MPI models. 

In contrast with angular, pT, and multiplicity dis- 
tributions of low Pt tracks traditionally used to test 
MPI models H, [l^l, we analyze events with high pT jets 
{pt > 15 GcV). Our approach complements the previous 
one since the MPI models have not been well tested in 
high Pt regimes, yet this kinematic region is the most 
important for searches for rare processes for which DP 
events are a potential background pTt - lsTf . 

This paper is organized as follows. In Sec. |lll we de- 
scribe the variables used in the analysis and motivate our 
choice of selection criteria. In Sec. IIIIl we describe the 



DO detector and the identification criteria for photons 
and jets. In Sec. llVi we describe the theoretical models 
used for comparison with data. In Sec. |Vl we discuss the 
corrections applied to the data in the cross section mea- 
surements and the related uncertainties. The measured 
cross sections and comparisons with some model predic- 
tions arc presented in Scc. lVIl In Sec. lVIIl we extract the 
fraction of DP events in the 7-I-2 jet final state as a func- 
tion of Pt*^- In Sec. IVIIIl we estimate the fractions of 
7-1-3 jet events occurring due to triple parton scattering, 
in bins of Pt*^^- Section HXl presents our conclusions. 

II. VARIABLES 

In this paper, we follow the notation used in our pre- 
vious analysis [2J] to distinguish between two classes of 
events. In events of the first class, the photon and all 
jets originate from the same single parton-parton inter- 
action (SP) with hard gluon bremsstrahlung in the initial 
or final state. In the second class, at least one of the jets 
originates from an additional parton interaction and thus 
we have at least two parton-parton interactions. 

To identify events with two independent parton-parton 
scatterings that produce a 7 -I- 3 jet final state, we use an 
angular distribution sensitive to the kinematics of the DP 
events [13l- We define an azimuthal angle between the 
Pt vector of the 7+lcading jet system and the pt sum of 
the two other jets: 

AS EE (P#, , (1) 
where = p^p + prf'^'^ and ~ Pt^^'^ + Pt^*^^ ■ 




FIG. 1: Diagram showing the pt vectors of the 7 + leading jet 
and jet2 jet3 systems in 7 + 3 jet events. 

Figure[l]shows the sum pT vectors of the 7+lcading jet 
and jet2 -I- jet3 systems in 7 + 3 jet events. In SP events, 
topologies with two radiated jets emitted close to the 
leading jet (recoiling against the photon direction in </>) 
are preferred and resulting in a peak at A5 = tt. How- 
ever, this peak is smeared by the effects of additional 
gluon radiation and detector resolution. For a simple 
model of DP events, with both the second and third jets 
originating from the second parton interaction, we have 
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exact pairwise balance in pT in both the 7 + jet and dijet 
system, and thus the A^* angle can have any value, i.e., 
we expect a uniform AS* distribution (35| . 

In this paper, we extend the study of DP interaction 
to the 7 + 2 jet events. In the presence of a DP inter- 
action the second jet in the event originates from a dijet 
system in the additional parton interaction and the third 
jet is either not reconstructed or below the pr threshold 
applied in the event selection. 

In the case of 7 + 2 jet events, we introduce a different 
angular variable, analogous to ([T|), to retain sensitivity to 
DP events. This variable is the azimuthal angle between 
the pt vector obtained by pairing the photon and the 
leading jet pT vectors and the second jet pT vector: 



A^ = A0 p^ 



.jct2 



(2) 



Figured] shows a diagram defining the pr of the two sys- 
tems in 7 -|- 2 jet events and the individual pt of the 
objects. The Ai/) distribution in 7 + 2 jet events has been 
used to estimate the DP fraction by the CDF Collabora- 
tion m. 




FIG. 2: Diagram showing the pr vectors of the 7 -I- leading jet 
system and p^°^'^ in 7 -I- 2 jet events. 



The Pt spectrum for jets from dijet events falls faster 
than that for jets, resulting from ISR and FSR in 7 + jet 
events, and thus the DP fractions should depend on the 
jet PT [il,[l-[ii,[2l. For this reason, the A</) dependent 
cross sections and the DP fractions in the 7 -I- 2 jet events 
are measured in three pi^*^ bins: 15 — 20, 20 — 25, and 
25 — 30 GeV. The AS* dependent cross section is mea- 
sured in 7 -|- 3 jet events (a subsample of the inclusive 7 
-I- 2 jet sample) in a single pi^*^ interval, 15 — 30 GeV. 
Such a measurement provides good sensitivity to the DP 
contribution, and discriminating power between different 
MPI models because the DP fraction in 7 -I- 3 jet events is 
expected to be higher than that in 7 + 2 jet events. This 
is expected since the second parton interaction will usu- 
ally produce a dijet final state, while the production of 
an additional jet in SP events via gluon bremsstrahlung 
is suppressed by the strong coupling constant Ug. 



III. DO DETECTOR AND DATA SAMPLES 

The DO detector is a general purpose detector de- 
scribed elsewhere in detail H^. Here we briefly describe 
the detector systems most relevant for this analysis. Pho- 
ton candidates are identified as isolated clusters of en- 
ergy deposits in the uranium and liquid-argon sampling 
calorimeter. The calorimeter consists of a central sec- 
tion with coverage in pseudorapidity |?7det| < 1-1 [13] and 
two end calorimeters covering up to |?7dct| ~ 4.2. The 
electromagnetic (EM) section of the calorimeter is seg- 
mented longitudinally into four layers, with transverse 
segmentation into cells of size Aiydet x A^dot = 0.1 x 0.1, 
except for the third layer, where it is 0.05 x 0.05. The 
hadronic portion of the calorimeter is located behind the 
EM section. The calorimeter surrounds a tracking sys- 
tem consisting of silicon microstrip and scintillating fiber 
trackers, both located within a solenoidal magnetic field 
of approximately 2 T. 

The events used in this analysis are required to pass 
triggers based on the identification of high Et clus- 
ters in the EM calorimeter with a shower shape con- 
sistent with that expected for photons. These triggers 
arc 100% efficient for photons with transverse momen- 
tum p7p > 35 GeV. To select photon candidates for our 
data sample, we use the following criteria [IJ, H^. EM 
objects are reconstructed using a simple cone algorithm 
with a cone size 7^ = 0.2 around a seed tower in 77 — 
space [3^. Regions with poor photon identification ca- 
pability and limited resolution (found at the bound- 
aries between calorimeter modules and between the cen- 
tral and end calorimeters) are excluded from the analy- 
sis. Each photon candidate is required to deposit more 
than 96% of its detected energy in the EM section of 
the calorimeter and to be isolated in the annular region 
between TZ = 0.2 and TZ = 0.4 around the center of 
the cluster: (i?f° - i?g°rc)/^j!Src < 0.07, where i?^^°t 
is the total (EM-|-hadronic) energy in the cone of ra- 
dius TZ = 0.4 and i?core the EM tower energy within 
a radius TZ = 0.2. The probability for candidate EM 
clusters to be spatially matched to a reconstructed track 
is required to be < 0.1%, where this probability is cal- 
culated using the spatial resolutions measured in data. 
We also require the energy-weighted EM cluster width 
in the finely-segmented third EM layer to be consistent 
with that expected for an electromagnetic shower. In ad- 
dition to calorimeter isolation, we apply track isolation, 
requiring that the scalar sum of the transverse momenta 
of tracks in an annulus of 0.05 < TZ < 0.4, calculated 
around the EM cluster direction, is less than 1.5 GeV. 

Jets arc reconstructed by clustering energy deposited 
in the calorimeter towers using the iterative midpoint 
cone algorithm with a cone size of 0.7. Jets must sat- 
isfy quality criteria that suppress background from lep- 
tons, photons, and detector noise effects. To reject back- 
ground from cosmic rays and W ^ ev decays, the missing 
transverse momentum, calculated as a vector sum of the 
transverse energies of all calorimeter cells, is required to 
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be less than 0.7 -p^. All pairs of objects (i, j) in the event 
(for example, photon and jet or jet and jet) are required 
to be separated by AR = ^(Ajyij)^ + {Acpi^y > 0.9. 

Each event must contain at least one photon in the 
pseudorapidity region |?7dot| < 1-0 or 1.5 < |?ydet| < 2.5 
and at least two (or three) jets with |r/dct| < 3.5. Events 
are selected with photon transverse momentum 50 < 
< 90 GeV, leading jet pt > 30 GeV, while the next- 
to-leading (second) jet must have pr > 15 GeV. If there 
is a third jet with pT > 15 GeV that passes the selec- 
tion criteria, the event is also considered for the 7 -I- 3 jet 
analysis. The higher p^ scale (i.e., the scale of the first 
parton interaction) , compared to the lower pT threshold 
required of the second (and eventual third) jet, results 
in a good separation between the first and second par- 
ton interactions of a DP event in momentum space. The 
reconstructed energy of each jet formed from calorime- 
ter energy depositions does not correspond to the actual 
energy of the jet particles which enter the calorimeter. 
It is therefore corrected for the energy response of the 
calorimeter, energy showering in and out the jet cone, 
and additional energy from event pile-up and multiple 
pp interactions. 

The sample of DP candidates is selected from events 
with a single reconstructed pp collision vertex. The col- 
lision vertex is required to have at least three associated 
tracks and to be within 60 cm of the center of the detector 
in the coordinate along the beam (z) axis. The probabil- 
ity for any two pp collisions occurring in the same bunch 
crossing for which a single vertex is reconstructed was 
estimated in and found to be < 10~^. 



IV. SINGLE AND MULTIPLE INTERACTION 
MODELS 

Monte Carlo (MC) samples are used for two purposes 
in this analysis. First, we use them to calculate recon- 
struction efficiencies and to unfold the data spectra to 
the particle level. Secondly, differential cross sections of 
7 -I- jet events simulated using different MPI models as 
implemented in the pythia and SHERPA [i^ event gen- 
erators are compared with the measured cross sections. 

There are two main categories of MPI models based 
on different sets of data used in the determination of 
the models parameters (it is customary to refer to dif- 
ferent models and to their settings as "tunes" ) . The two 
categories, "old" and "new," correspond to different ap- 
proaches in the treatment of MPI, ISR and FSR, and 
other effects [l^ [l^l . The main difference between the 
"new" [ll| and the "old" models is the implementation 
of the interplay between MPI and ISR, i.e., considering 
these two effects in parallel, in a common sequence of 
decreasing pT values. In the "old" models, ISR and FSR 
were included only for the hardest interaction, and this 
was done before any additional interactions were consid- 
ered. In the "new" models, all parton interactions in- 
clude ISR and FSR separately for each interaction. The 



new models, especially those corresponding to the Peru- 
gia family of tunes jlOj , also allow for a much wider set of 
physics processes to occur in the additional interactions. 
A detailed description of the different pythia MPI mod- 
els can be found elsewhere [l^, [HI . Here we provide a 
brief description of the models considered in our anal- 
ysis. They include: Perugia-0 (PO, the default model 
in the Perugia family [13); P-hard and P-soft, which 
explore the dependence on the strength of ISR/FSR ef- 
fects, while maintaining a roughly-consistent MPI model 
as implemented in the PO tune; P-nocr, which excludes 
any color reconnections in the final state; P-X and P-6, 
which are PO modifications based on the MRST LO* and 
CTEQ6L1 PDF sets, respectively (PO uses CTEQ5L as 
a default). We also compare data with predictions de- 
termined using tunes A and DW as representative of the 
"old" MPI models. 

The measured cross sections are also compared with 
predictions obtained from the SHERPA event generator, 
which also contains a simulation of MPI. Its initial mod- 
elling was similar to tunc A from pythia 0, but it has 
evolved and now is characterized, in particular, by (a) 
showering effects in the second interaction and (b) a com- 
bination of the CKKW merging approach with the MPI 
modeling [13, SH . Another distinctive feature of SHERPA 
is the modeling of the parton-to-photon fragmentation 
contributions through the incorporation of QED effects 
into the parton shower p^ . 

The data are also compared with models without MPI, 
in which the photon and all the jets are produced exclu- 
sively in SP scattering. Such events are simulated in 
both PYTHIA and SHERPA. In PYTHIA, only 2 2 di- 
agrams are simulated, resulting in the production of a 
photon and a leading jet. With the MPI event generation 
switched off, all the additional (to the leading) jets are 
produced in the parton shower development in the initial 
and final states. We refer to such SP events as "pythia 
SP" events. In SHERPA, up to two extra partons (and 
thus jets) are allowed at the matrix clement level in the 
2 — > {2,3,4} scattering, but jets can also be produced 
in parton showers. To provide a matching between the 
matrix-element partons and parton shower jets, we follow 
the recommendation provided in [Z^] and choose the fol- 
lowing "matching" parameters: the energy scale Qo = 30 
GeV and the spatial scale D = 0.4, where D is taken to 
be of the size of the photon isolation cone [1^. This is 
the default scheme for the production of 7 -I- jet events 
with and without MPI simulation. The set of events 
produced without MPI simulation within this scheme is 
called "SHERPA- 1 SP". We study the dependence of the 
measured DP fractions on the scale choice in the SHERPA 
SP models in Sec. IVIII by setting the matching scale 
Qo equal to 20 and to 40 GeV (sets "sherpa-2" and 
"SHERPA-3" respectively). For completeness we consider 
SHERPA SP events in which all of the extra jets are pro- 
duced (as in pythia) in the parton shower with only a 
2 — 2 matrix element, and caff this set "sherpa-4 SP". 
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V. DATA ANALYSIS AND CORRECTIONS 
A. Background studies 

The main background to isolated photons comes from 
jets in which a large fraction of the transverse momen- 
tum is carried by photons from tt", 77, or decays. The 
photon-enhancing criteria described in Sec. IIIII are de- 
veloped to suppress this background. The normalized 
AS and Acj) distributions are not very sensitive to the 
exact amount of background from events with misiden- 
tified photons. To estimate the photon fractions in the 
Acj) bins, we use the output of two neural networks (NN) 
[38j . These NNs are constructed using the jetnet pack- 
age [4^ 1 and are trained to discriminate between photon 
and EM-jets in the central and end calorimeter regions 
using calorimeter shower shape and track isolation vari- 
ables [11] . The distribution of the photon NN output for 
the simulated photon signal and for the dijet background 
samples are fitted to data in each A(f) bin using a max- 
imum likelihood fit [45| to obtain the fractions of signal 
events in the data. To obtain a more statistically signifi- 
cant estiniate of the photon purity in the Acj) bins, we use 
a single p^^*^ bin: 15 < p!^^^ < 30 GeV. The fit resuhs 
show that the 7-l-jet signal fractions in all A(j> bins agree 
well with a constant value, 0.69 ± 0.03 in the central and 
0.71 ± 0.02 in the end calorimeter regions. 

The sensitivity of the AS and Acj) distributions to this 
background from jets is also examined by considering two 
data samples in addition to the sample with the default 
photon selections: one with relaxed and another with 
tighter track and calorimeter isolation requirements. Ac- 
cording to MC estimates, in those two samples the frac- 
tion of background events should either increase or de- 
crease by (30 — 35)% with respect to the default sam- 
ple. We study the variation of the AS and A0 normal- 
ized cross sections in data by comparing the relaxed and 
tighter data sets with the default set. We find that the 
cross section variations are within 5%. 



B. Efficiency and unfolding corrections 

To select j + 2 jet and 7-I-3 jet events, we apply the se- 
lection criteria described in Sec. IIIII The selected events 
are then corrected for selection efficiency, acceptance, and 
event migration effects in bins of AS and Acj). These cor- 
rections are calculated using MC events generated using 
PYTHIA with tune PO, as discussed in Sec. |lVl The 
generated MC events are processed through a GEANT- 
based [i^ simulation of the DO detector response. These 
MC events are then processed using the same reconstruc- 
tion code as used for the data. We also apply additional 
smearing to the reconstructed photon and jet pr so that 
the resolutions in MC match those observed in data. The 
reconstructed AS and A(p distributions in the simulated 
events using the PO tune are found to describe the data. 
In addition to the simulation with the default tune PO, 



we have also considered PO MC events that have been 
reweighted to reproduce the pr distributions in data. 
After such reweighting, the reconstructed AS* and A<j> 
distributions give an excellent description of the data. 
Figures [3] and |4] show the normalized distributions as a 
function of AS* for 7 -I- 3 jet and A(f) for 7 + 2 jet events 
(for the p-!^*^ bin 15 — 20 GeV, chosen as an example) in 
data and in the reweighted MC. 
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FIG. 3: Normalized AS distribution for data and for the 
reweighted MC sample in the range 15 < p^*^ < 30 GeV. 
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FIG. 4: Normalized A(j} distribution for data and for the 
reweighted MG sample in the range 15 < Py*^ < 20 GeV. 

Three sets of corrections are applied in data to ob- 
tain the differential cross sections which we then com- 
pare with the various MPI models. We apply them to 
correct for detector and reconstruction inefficiencies and 
for bin migration effects. The first correction deals with 
the possibility that, due to the detector and reconstruc- 
tion effects, our selected event sample may contain events 
which would fail the selection criteria at the particle level. 
The data distributions are corrected, on a bin-by-bin ba- 
sis, for the fraction of events of this type. We also apply a 
correction for events which fail the selection requirement 
at the reconstruction level. Systematic uncertainties are 
assigned on these two correction factors to account for 
uncertainties on the photon and jet identification, the jet 
energy scale (JES) and resolution, and vary in the AS 
(Acj)) bins up to 12% (18%) in total. They are domi- 
nated by the JES uncertainties. The overall corrections 
obtained with the default PO and with the reweighted 
MC samples agree within about 5% for most AS" and A(j> 
bins and differ by at most 25%. Since we are measur- 
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ing normalized cross sections, the absolute values of the 
corrections are not important, and we need only their 
relative dependence on AS and A(j). 

The third correction accounts for the migration of 
events between different bins of the AS* and distribu- 
tions, which is caused by the finite photon and jet angu- 
lar resolutions and by energy resolution effects, and can 
change the pt ordering of jets between the reconstruc- 
tion and the particle level. To obtain the AS and A0 
distributions at the particle level, we follow the unfold- 
ing procedure described in the Appendix, based on the 
Tikhonov regularization method j47H5C)l |. The bin sizes 
for the AS* and A0 distributions are chosen to have sensi- 
tivity to different MPI models (which is largest for small 
AS* and A(p angles) while keeping good statistics and the 
bin-to-bin migration small. The statistical uncertainties 
(<5stat) arc in the range (10 — 18)%. They arc due to the 
procedure of regularized unfolding and take into account 
the correlations between the bins. The correlation factor 
for adjacent bins in the unfolded distributions is about 
(30 — 45)%, and it is reduced to w 10% for other (next- 
to-adjacent) bins. To validate the unfolding procedure, a 
MC closure test is performed. Wc compare the unfolded 
MC distribution to the true MC distribution and find 
that they agree within statistical uncertainties. 



VI. DIFFERENTIAL CROSS SECTIONS AND 
COMPARISON WITH MODELS 

In this section, we present the four measure- 
ments of normalized differential cross sections, 
{l/a~f3j)daj3j /dAS in a single p^^*^ bin (15 — 30 
GeV) for 7 -I- 3 jet events and {l/ay2j)d(Ty2j /dA(j) in 
three p^^*^ bins (15 - 20, 20 - 25, and 25 - 30 GeV) for 
7 -|- 2 jet events. The results arc presented numerically 
in Tables H] - IIVI as a function of AS* and A<j>, with the 
bin centers estimated using the theoretical predictions 
obtained using the PO tune. The column 'Wdata" shows 
the number of selected data events in each AS* (Ac/)) bin 
at the reconstruction level. The differential distributions 
decrease by about two orders of magnitude when moving 
from the AS (Acj)) bin 2.85 — 3.14 radians to the bin 
0.0 — 0.7 radians and have a total uncertainty ((5tot) 
between 7% and 30%. Here (Stot is defined as a sum 
in quadrature of statistical ((5stat) and systematical 
((5syst) uncertainties. It is dominated by systematic 
uncertainties. The sources of systematic uncertainties 
are the JES (2 — 17%), largest at the small angles, 
unfolding (5 — 18%), jet energy resolution simulation in 
MC events (1 — 7%), and background contribution (up 
to 5%). 

The results are compared in Figs. [SHHl to predictions 
from different MPI models implemented in pythia and 
SHERPA, as discussed in Sec. |lVl We also show predic- 
tions of SP models in pythia and SHERPA (sherpa-1 
model). In the QCD NLO predictions, only final states 
with a 7 + 2 jet topology arc considered, and thus for 



the direct photon production diagrams, we should have 
Acj) = TT. The Acj) angle may differ from tt due to pho- 
ton production through a parton-to-photon fragmenta- 
tion mechanism. Even if we take into account this pro- 
duction mechanism, which is included in the jetphox 
[Slj NLO QCD calculations, only the two highest A(j) 
bins receive significant contributions. 

Figs. [5]-[5] show the sensitivity of the two angular vari- 
ables AS and A^ to the various MPI models, with predic- 
tions varying significantly and differing from each other 
by up to a factor 2.5 at small AS and Ac/), in the re- 
gion where the relative DP contribution is expected to 
be the highest. The sensitivity is reduced by the choice 
of SP model, for which we derive an upper value of 25% 
comparing the ratios of predictions from various models 
(pythia, SHERPA-2, -3, -4). This upper value is con- 
siderably smaller than the difference between the various 
MPI models. 

Tables IVl and I VII are complementary to Figs. [5] -[5] and 
show the x^/ndf values of the agreement between theory 
and data for each model. Here ndf stands for the number 
of degrees of freedom (taken as the number of bins, iVbins 
minus 1), and is calculated as 

,3, 

i—l i.unc 

where Di and represent the cross section values in 
the i-th bin of data and a theoretical model respectively, 
while (Sf ujjc is the total uncorrelated uncertainty in this 
bin. The latter is composed of the uncertainties for the 
corrections in the unfolding procedure (Sec. IVSp . the 
statistical uncertainties of the data (5stat and the theoret- 
ical model. The uncorrelated uncertainty is always 
larger than all remaining correlated systematic uncer- 
tainties. Since small angles {AS{A(j>) < 2) are the most 
sensitive to DP contributions, we calculate the /ndf 
separately for these bins. From Figs. [5]- 18] and Tables IVl 
and lVIi we conclude: (a) the predictions derived from SP 
models do not describe the measurements; (b) the data 
favor the predictions of the new pythia MPI models (PO, 
P-hard, P-6, P-X, P-nocr) and to a lesser extent SO and 
SHERPA with MPI; and (c) the predictions from tune A 
and DW MPI models are disfavored. 



VII. FRACTIONS OF DOUBLE PARTON 
EVENTS IN THE 7 + 2 JET FINAL STATE 

The comparison of the measured cross section with 
models (Sec. IVip shows clear evidence for DP scatter- 
ing. We use the measurement of the differential cross 
section with respect to A0 and predictions for the SP 
contributions to this cross sections in different models to 
determine the fraction of 7 -I- 2 jet events which originate 
from DP interactions as a function of the second par- 
ton interaction scale (p^^'^) and of A0. Due to ISR and 
FSR effects the pT balance vectors of each system may be 
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TABLE I: Measured normalized differential cross sections 
(l/crT,3j )dcr^:y/dAS for 15 < pj^" < 30 GeV. 



AS bin i^S) Ndata Normalized Uncertainties (%) 



(rad) 


(rad) 




cross section Sstat 




Stot 


0.00 


- 0.70 


0.36 


495 


2.97 X 10" 


■2 


11.3 


14.7 


18.6 


0.70 


- 1.20 


0.97 


505 


4.74 X 10" 


2 


12.3 


15.6 


19.9 


1.20 


~ 1.60 


1.42 


498 


5.80 X 10" 


■2 


13.4 


15.8 


20.7 


1.60 


- 2.15 


1.90 


1315 


1.11 X 10" 


-1 


7.5 


15.3 


17.0 


2.15 


- 2.45 


2.32 


1651 


2.38 X 10" 


■1 


6.0 


12.0 


13.4 


2.45 


- 2.65 


2.56 


1890 


4.04 X 10" 


-1 


5.6 


13.6 


14.7 


2.65 


- 2.85 


2.76 


3995 


8.59 X 10" 


■1 


3.2 


5.6 


6.4 


2.85 


- 3.14 


3.02 


12431 


1.89 X 10" 




1.0 


13.0 


13.0 



TABLE III: Measured normalized differential cross section 
(l/(T^2j )dcr^2j/dA<^ for 20 < p^*"^ < 25 GeV. 



Acj) bin 


{A(j>) iVdata 


Normalized 


Uncertainties (%) 


(rad) 


(rad) 




cross section 5stat 




Stot 


0.00 - 0.70 


0.35 


388 


1.17 X 10"^ 


12.5 


23.2 


26.4 


0.70 - 1.20 


0.96 


358 


1.75 X 10"^ 


17.7 


22.2 


28.5 


1.20 - 1.60 


1.42 


489 


3.29 X 10"^ 


15.6 


17.0 


23.1 


1.60 - 2.15 


1.92 


1848 


9.84 X 10"^ 


6.2 


13.8 


15.1 


2.15 - 2.45 


2.33 


2682 


2.80 X 10"^ 


4.6 


8.2 


9.4 


2.45 - 2.65 


2.56 


3208 


5.21 X 10"^ 


4.5 


7.1 


8.4 


2.65 - 2.85 


2.77 


5404 


9.01 X 10"^ 


3.2 


7.3 


8.0 


2.85 - 3.14 


3.02 


15901 


1.88 X 10° 


1.0 


10.8 


10.8 



TABLE IL Measured normalized differential cross sections TABLE IV: Measured normalized differential cross section 
{l/a^2j)da^2j /dA(j) for 15 < p^^" < 20 GeV. (l/cr^2j)rfcr^2j/rfA(^ for 25 < p^^*^ < 30 GeV. 



A(j> bin 


{Acj>) A^data 


Normalized 


Uncertainties (%) 


A<j!> bin 


(A0) 


A^data 


Normalized 


Uncertainties (%) 


(rad) 


(rad) 




cross section (5stat 




Stot 


(rad) 


(rad) 




cross section Ssta.t 


(5syst 


Stot 


0.00 - 0.70 


0.36 


1028 


2.49 X 10"^ 


9.4 


19.1 


21.3 


0.00 - 0.70 


0.32 


158 


6.82 X 10"^ 


16.1 


19.8 


25.5 


0.70 - 1.20 


0.96 


822 


3.06 X 10-2 


11.8 


20.3 


23.4 


0.70 - 1.20 


0.94 


155 


1.11 X 10"2 


20.9 


16.4 


26.6 


1.20 - 1.60 


1.42 


1149 


5.68 X 10"2 


9.6 


15.5 


18.2 


1.20 - 1.60 


1.45 


190 


1.87 X 10"2 


24.0 


17.9 


30.0 


1.60 - 2.15 


1.92 


3402 


1.29 X 10"^ 


4.9 


11.5 


12.5 


1.60 - 2.15 


1.92 


910 


7.00 X 10-2 


7.0 


15.9 


17.4 


2.15 - 2.45 


2.32 


4187 


3.06 X 10"^ 


4.5 


9.5 


10.5 


2.15 - 2.45 


2.32 


1683 


2.50 X 10"^ 


5.0 


8.6 


9.9 


2.45 - 2.65 


2.56 


5239 


5.88 X 10"^ 


4.0 


6.3 


7.4 


2.45 - 2.65 


2.57 


2155 


4.93 X 10"^ 


4.5 


8.9 


10.0 


2.65 - 2.85 


2.76 


8246 


9.43 X 10"^ 


3.0 


6.8 


7.5 


2.65 - 2.85 


2.77 


3894 


9.09 X 10-1 


3.1 


7.5 


8.1 


2.85 - 3.14 


3.01 


20337 


1.63 X 10° 


1.1 


12.3 


12.3 


2.85 - 3.14 


3.03 


12332 


2.01 X 10° 


1.0 


10.2 


10.2 



non-zero and have an arbitrary orientation with respect 
to each other [sH , which leads to a uniform Ac/) distribu- 
tion for DP events. 

Using the uniform distribution as the DP model tem- 
plate and the SHERPA-1 prediction as the SP model tem- 
plate, we can fit the A(f> distributions measured in data 
and obtain the fraction of DP events from a maximum 
Ukelihood fit [i^. We repeat this procedure in three in- 
dependent ranges of p^^'^. The distributions in data, SP, 
and DP models, as well as a sum of the SP and DP dis- 
tributions, weighted with their respective fractions, arc 
shown in Figs. [9] - [TT] for the three p^^'^ intervals. The 
sum of the SP and DP predictions reproduces the data 
well. The measured DP fractions ifjp'') are presented in 
Table IVnl 

The uncertainties in the DP fractions are due to the 
fit, the total (statistical plus systematic) uncertainties on 
the data points, and the choice of SP model. The effect 
from the second source is estimated by varying all the 
data points simultaneously up and down by the total ex- 
perimental uncertainty (Stot)- The uncertainties due to 
the SP model are estimated by considering SP models 
(Sec. IIVP that arc different from the default choice of 
SHERPA-1: -2, -3, -4, as well as pythia SP predictions. 
The measured DP fractions with all sources of uncertain- 
ties in each p^*^ bin are summarized in Table IVIII The 

DP fractions in 7-I-2 jet events, fj^'', decrease as a func- 
tion of pj^*^ from (11.6 ± 1.0)% in the bin 15 - 20 GeV, 



to (5.0 ± 1.2)% in the bin 20 - 25 GeV, and (2.2 ± 0.8)% 
in the bin 25 — 30 GeV. The estimated DP fraction in 
7 -I- 2 jet events selected with pj > 16 GeV and p)^' > 8 
GeV from the CDF Collaboration [H is 14+^%, which 
is in qualitative agreement with an extrapolation of our 
measured DP fractions to lower jet pt- 

The DP fractions shown in Table IVIII are integrated 
over the entire region < A(j) < tt. However, from Figs. [5] 
- [Til the fraction of DP events is expected to be higher 
at smaller Ac/) . To determine the fractions as a function 
of A(j) , we perform a fit in the different A^ regions by 
excluding the bins at high Acf) ; specifically, by consider- 
ing the Act) regions - 2.85, - 2.65, - 2.45, - 2.15, 
and — 1.60. The DP fractions for these A0 regions are 
shown in Table IVIIII for the three pJ^*^ intervals. The 
DP fractions with total uncertainties as functions of the 
upper limit on A(j> (A(/)max) for all the pj^'^ bins arc also 
shown in Fig. 1121 As expected, they grow significantly 
towards the smaller angles and are higher for smaller p^^*^ 
bins. 



VIII. FRACTIONS OF TRIPLE PARTON 
EVENTS IN THE 7 + 3 JET FINAL STATE 

In this section, we estimate the fraction of 7 + 3 jet 
events from triple parton interactions (TP) in data as a 
function of p^^*^. In 7 -|- 3 jet TP events, the three jets 
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TABLE V: The results of a test of the agreement between data points and theory predictions for the AS (7 + 3 jet) and 

A(j!> (7 + 2 jet) distributions for 0.0 < AS(A<j!)) < vr rad. Values are x^/ndf. 

Variable p^^*^ SP model MPI model 

(GeV) PYTHIA SHERPA A DW SO PO P-nocr P-soft P-hard P-6 P-X sherpa 

A5 15 - 30 7?7 6i) 15.6 21.4 2.2 0.4 05 2^9 05 Ol 05 

A(^ 15 - 20 16.6 11.7 19.6 27.7 1.6 0.5 0.9 1.6 0.9 0.6 0.8 1.2 

A(^ 20 - 25 10.2 5.9 4.0 7.9 1.1 0.9 1.4 2.1 1.1 1.3 1.5 0.4 

A(^ 25 - 30 7.2 3.5 2.8 3.0 2.4 1.1 1.1 3.7 0.2 1.3 1.9 0.7 



TABLE VI: The results of a test of the agreement between data points and theory predictions for the AS (7 + 3 jet) and 
the A(/> (7 + 2 jet) distributions for AS(A(^) < 2.15 rad. Values are /ndf . 



Variable 


ict2 

Ft 


SP model 










MPI model 












(GeV) 


PYTHIA 


SHERPA 


A 


DW 


SO 


PO P 


-nocr 


P-soft P-hard 


P-6 


P-X 


SHERPA 


AS 


15-30 


10.9 


11.3 


31.0 


42.9 


3.4 


0.4 


0.5 


4.9 


0.9 


0.5 


0.4 


2.6 


A<^ 


15-20 


30.2 


26.0 


40.7 


61.1 


2.2 


0.9 


1.6 


1.5 


1.2 


1.2 


1.0 


2.4 


A(^ 


20 - 25 


15.4 


12.1 


6.8 


18.0 


1.0 


1.8 


2.7 


1.7 


1.5 


2.5 


2.4 


0.6 


A(^ 


25-30 


7.1 


5.3 


1.3 


5.6 


1.6 


1.1 


1.0 


2.1 


0.3 


1.4 


1.6 


0.5 



CO 

<l 

T3 



: D0, L = 1 .0 fb^' 



10': 



o 

(U 



(0 
Q 



Data 

PYTHIA, tune A 
PYTHIA, tune DW 
PYTHIA, tune SO 
PYTHIA, tune PO 
SHERPA, with MPI 
PYTHIA, no MPI 
SHERPA, no MPI 
j Total uncertainty 



50 < p^^ < 90 GeV 
p!^'" > 30 GeV 
15<p!^°'^<30 GeV 
p';''>15GeV 



I 



□ 



_l_ 



1.4 r 

ir A . ^ ^ ^ A.:* 

0.8- ^. ....f.... » I I' 

0.6 r ♦ * * 

, * I , , , , I , , , I , , , , I , , , , I , , , 

0.5 1 1.5 2 2.5 

AS 



3 

(rad) 



-e- 
<l 

■D 



T3 



10 ^ 



10 



lO' 



00, L =1.0fb"' 

int 

• Data 

i PYTHIA, tune A 
PYTHIA, tune DW 
PYTHIA, tune SO 

A PYTHIA, tune PO 

* SHERPA, with IVIPI 
o PYTHIA, no MPI 

□ SHERPA, no MPI 
: Total uncertainty 

I I I 



50 < p^ < 90 GeV 
30 GeV 

20 GeV 



P t 



15 <p'^' 




2.5 3 

A(|) (rad) 



FIG. 5: Normalized differential cross section in 7 + 3 jet 
events, {l/a~,3j)da-,3j /dAS, in data compared to MC models 
and the ratio of data over theory, only for models including 
MPI, in the range 15 < p^^^ < 30 GeV. 



FIG. 6: Normalized differential cross section in 7 + 2 jet 
events, {l/a-,2j)da~,2j /dA<j}, in data compared to MC models 
and the ratio of data over theory, only for models including 
MPI, in the range 15 < pj^*^ < 20 GeV. 



TABLE VII: Fractions of DP events (%) with total uncertain- 
ties for < A(p < TT in the three p^*^ bins. 



Pt 

(GeV) 


(GeV) 


•/dp 
(%) 


Uncertainties (in %) 
Fit Stot SP model 


15-20 


17.6 


11.6 ± 1.4 


5.2 8.3 6.7 


20 - 25 


22.3 


5.0 + 1.2 


4.0 20.3 11.0 


25-30 


27.3 


2.2 + 0.8 


27.8 21.0 17.9 



come from three different parton interactions, one 7+ jet 
and two dijet final states. In each of the two dijet events, 
one of the jets is either not reconstructed or below the 
15 GeV Pt selection threshold. 

In our previous study of DP 7 + 3 jet events [l^, we 
built a data-driven model of inclusive DP interactions 
(MixDP) by combining 7+jet and dijet events from data, 
and obtaining the 7 + 3 jet + X final state. However, 
since each component of the MixDP model may contain 
two (or more) jets, where one jet is due to an additional 
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TABLE VIII: DP fractions (%) in data as a function of the A(j) interval for three Py*^ bins. 
Py''^ interval (rad) 

(GeV) O-TT - 2.85 0- 2.65 0- 2.45 - 2.15 0- 1.6 
15-20 11.6 ±1.4 18.2 ± 2.4 25.0 ± 2.9 33.7 ± 3.8 45.0 ± 5.5 47.4 ±11.4 
20 - 25 5.0 ±1.2 9.4 ±1.2 13.4 ±2.1 19.6 ±3.1 28.1 ±4.3 63.7 ±17.2 
25 - 30 2.2 ±0.8 3.8 ± 1.3 5.0 ±1.5 6.2 ± 2.2 9.8 ± 4.5 27.8 ±11.5 



D0, L, = 1.0 fb"' 

int 

• Data 

i PYTHIA, tuneA 

PYTHIA, tune DW 
■ PYTHIA, tune SO 
A PYTHIA, tune PO 
O SHERPA, with MPI 
O PYTHIA, no MPI 
□ SHERPA, no MPI 
; Total uncertainty 



50 < < 90 GeV 

p'°" > 30 GeV 

20 < p'°'^ < 25 GeV 




0.5 



1.5 



2.5 3 

A(|) (rad) 



-e- 
<l 



T3 



Data 

PYTHIA, tuneA 
PYTHIA, tune DW 
PYTHIA, tune SO 
PYTHIA, tune PO 
SHERPA, with MPI 
PYTHIA, no MPI 
SHERPA, no MPI 
; Total uncertainty 



50 < p^ < 90 GeV 



|iat1 , 

25 <p';" 



c 30 GeV 




FIG. 7: Normalized differential cross section in 7 ± 2 jet 
events, {l/ary2j)da-,2j /dA(t>, in data compared to MC models 
and the ratio of data over theory, only for models including 
MPI, in the range 20 < p^^*^ < 25 GeV 



FIG. 8: Normalized differential cross section in 7 ± 2 jet 
events, {l/(j-y2j)dG^,2i /dAcj), in data compared to MC models 
and the ratio of data over theory, only for models including 
MPI, in the range 25 < p^^'^ < 30 GeV. 



parton interaction, the model simulates the properties of 
"double plus triple" parton interactions. Therefore, the 
"DP" fractions found earlier in the 7 + 8 jet data (shown 
in Table III of [24jj ) take into account a contribution from 
TP interactions as well. These fractions are also shown in 
the second column of Table HXl Thus, if we calculate the 
fractions of TP events in the MixDP sample, defined as 
/tp'"'^'^' '^^'^ calculate the TP fractions in the 7±3 jet 
data, f'^p ■' , as 

f73j _ /■dp+tp r73j rA\ 
Jtp ~ Jtp ■ ./dp+tp' 

where /^p+tp ^^'^ fraction of DP+TP events in the 
7 ± 3 jet sample. Figure [13] shows two possible ways in 
which a DP event and an SP event can be combined to 
form a 7 + 3 jet event which is a part of the MixDP 
sample, with details on the origin of the various parts of 
the event given in the caption. Contributions from other 
possible MixDP configurations are negligible (< 1%). In 
[24| . we calculated how often each component. Type I and 
II, is found in the model. Table HXl shows that the events 




0.5 1 1.5 2 2.5 3 

Alt) (rad) 



FIG. 9: A(j!> distribution in data, SP, and DP models, and 
the sum of the SP and DP contributions weighted with their 
fractions for 15 < p^^" < 20 GeV. 
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FIG. 10: Same as in Fig. [H but for 20 < pi^'^ < 25 GeV. 



Type I 



Type II 





FIG. 13: Two possible combinations of events present in the 
MixDP sample that actually represent the contribution from 
triple scattering events in the 7 + 3 jet final state: (a) a 7 + 
1 jet event mixed with a double-dijet DP event, where one jet 
of each dijet is lost or not reconstructed (Type I); (b) a DP 
event in the (7+ 1 jet)+dijet final state with one jet from the 
dijet lost, mixed with a dijet event with one jet lost (Type 
II). Dashed lines correspond to the lost jets. 



of Type II arc dominant in all bins. Thus, the fraction of 
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FIG. 11: Same as in FigM but for 25 < p't < 30 GeV. 
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FIG. 12: Fractions of DP events with total uncertainties in 
7 + 2 jet final state as a function of the upper limit on for 
the three Py'^ intervals. 



TABLE IX: Fractions of DP+TP events with total uncer- 
tainties in 7 + 3 jet data (/^p+tp) ^^'^ fractions of Type I (II) 



events in the data-driven DP model {3 



p^^^ bins. 



Typo I(II) 



) in the three 



Ft 

(GeV) 


Jdp-l-tp 
(%) 


-FTypc I 


-Frypo II 


15-20 


46.6 ±4.1 


0.26 


0.73 


20 - 25 


33.4 ±2.3 


0.22 


0.78 


25-30 


23.5 ±2.7 


0.14 


0.86 



TP configurations (Fig.HS]) in the MixDP model, ftp^^"^, 
can be calculated as 



/.dp-l-tp 

Jtp 



''Typo II ■ /dp 



Typo I ■ J(jp 



(5) 



where fj^^ and /j^ are the fractions of events with DP 
scattering resulting in 7±2 jet and dijet final states. We 
separately analyze each of the event types of Fig. [T31 The 
fraction of events having a second parton interaction with 
a dijet final state with cross section cr-'-' can be defined 
using the effective cross section (Toff as f^^ = a^^ / (2crcff ). 
The cross section for a DP scattering producing two di- 



jet final states can be presented then as o'dp"'-' 



dp 



0, ij. The fraction /j^ is estimated using dijet events 
simulated with pythia. We calculate the jet cross sec- 
tions cr^^ for producing at least one jet in the three pr 
bins with \rf'^^\ < 3.5. The effective cross section acs is 
taken as an average of the CDF [2l| and DO [l^l mea- 
surements, (Tgg°'^^^'^ = 15.5 mb. The determined fractions 
are shown in the third column of Table [Xj We assume 
that the estimates, done at the particle level, are also ap- 
proximately correct at the reconstruction level. We take 
an uncertainty on these numbers Sf^^ = f^^. 

The fractions of the 7+2 jet events in which the second 
jet is due to an additional parton scattering arc estimated 
in the previous section and arc much higher than /^p. 
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TABLE X: Fractions of DP events in 7 + 2 jet {f2p^) and dijet 
(/dp) final states as well the fraction of TP configurations in 



« 0.08 



\.ldp 

the MixDP model (/tp''+*'') 



in the three p^*^ bins 



(GeV) 


f7^J 
Jdp 

(%) 


Jdp 
(%) 


/■dp+tp 
Jtp 

(%) 


15-20 


15.9 ±2.2 


0.50 ±0.50 


11.7 ± 1.9 


20 - 25 


7.8 ±2.0 


0.17 ±0.17 


6.1 ± 1.8 


25-30 


4.2 ± 1.3 


0.07 ±0.07 


3.6 ± 1.2 



However, since we estimate the TP fraction in data at the 
reconstruction level, we repeat the same fitting procedure 
used for the extraction of /^p"' from the Acj) distributions 
in the reconstructed data and SP 7 + 2 jet MC events. 
The results of the fit in the three p^^*^ intervals are sum- 
marized in the second column of Table [X] Here the total 
uncertainties (5tot are due to the statistical and system- 
atic uncertainties shown in Tables ITll- II V[ but excluding 
the uncertainties from the unfolding. 

By substituting f^-'^ and fjp'' into Eq. (O , we calculate 

the TP fractions ftp^^^ in the MixDP model. They arc 
shown in the last column of Table [Xj The TP fraction in 
the similar data-driven MixDP model in the CDF anal- 
ysis (2]| for (JES uncorr ected) 5 < ^^'^ < 7 GeV was 
estimated as 17lg%, i.e., a value that is higher, on av- 
erage, than our TP fractions measured at higher jet px, 
but in agreement with an extrapolation of our observed 
trend to lower jet pT- 

By substituting /t'jf'''''' and the DP+TP fractions in 

7±3 jet data fj^+tp from ^ into Eq. (gD, we get the TP 

fractions in the 7 ± 3 jet data, f^p'', which are shown in 
the second column of Table IXII They are also presented 
in Fig. [131 The pure DP fractions, /^p"' , can then be 

obtained by subtracting the TP fractions f^p-' from the 
inclusive DP+TP fractions fjpl^p- 

The last column of Table IXII shows the ratios of the 
TP to DP fractions /tjp^V/dp"' in 7 + 3 jet events. Since 
the probability of producing each additional parton scat- 
tering with a dijet final state is expected to be directly 
proportional to cr^^/cToff, the ftp^ / f2'p ratio should be 
approximately proportional to the jet cross section a^^ , 
and drop correspondingly as a function of the jet pT- 
This trend is confirmed in Table IXll 



TABLE XI: Fractions of TP events (%) and the ratio of 
TP/DP fractions in the three p^*^ bins of 7 + 3 jet events. 
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FIG. 14: Fractions of TP events with total uncertainties in 
7 + 3 jet final state as a function of p^*'^. 



IX. SUMMARY 

We have studied the azimuthal correlations in 7 + 
3 jet and 7 + 2 jot events and measured the normal- 
ized differential cross sections (l/cr-y3j)(icr-y3j/(iAS' and 
{1 / aj2j)daj2j / dA(j) in three bins of the second jet pT- 
The results are compared to different MPI models and 
demonstrate that the predictions of the SP models do 
not describe the measurements and an additional contri- 
bution from DP events is required to describe the data. 
The data favor the predictions of the new pythia MPI 
models with p^-ordcrcd showers, implemented in the Pe- 
rugia and SO tunes, and also SHERPA with its default 
MPI model, while predictions from previous pythia MPI 
models, with tunes A and DW, are disfavored. 

We have also estimated the fractions of DP events in 
the 7 + 2 jet samples and found that they decrease in the 
bins of pj^*^ as (11.6 + 1.0)% for 15-20 GeV, (5.0 + 1.2)% 
for 20-25 GeV, and (2.2+0.8)% for 25-30 GeV. Finally, 
for the first time, we have estimated the fractions of TP 
events in the 7 + 3 jet data. They vary in the Pj.^^ bins 
as (5.5 + 1.1)% for 15 - 20 GeV, (2.1 + 0.6)% for 20 - 25 
GeV, and (0.9 + 0.3)% for 25 - 30 GeV. 

The measurements presented in this paper can be used 
to improve the MPI models and reduce the existing theo- 
retical ambiguities. This is especially important for stud- 
ies in which a dependence on MPI models is a significant 
uncertainty (such as the top quark mass measurement), 
and in searches for rare processes, for which DP events 
can be a sizable background. 
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procedure may introduce a bias (6). We optimize the 
regularization by finding a balance between Jgtat and b 
according to the following criterion: we minimize the fol- 
lowing function of i5stat and b in the first two bins, — 0.7 
and 0.7-1.2: 



X. APPENDIX 

In this appendix we discuss the unfolding procedure 
used to correct the measured AS and A(j) distributions 
to the particle level in order to obtain a nonparamet- 
ric estimate of the true AS* and A<^ distributions from 
the measured (reconstructed) distribution taking into ac- 
count possible biases and statistical uncertainties. We 
use the following approach to extract the desired distri- 
butions from our measurements. The observed distribu- 
tion is the result of the convolution of a resolution func- 
tion with the desired distribution at the particle level. 
After the discretization in AS/A<f> bins, the resolution 
function is a smearing matrix, and distributions on both 
particle level and reconstruction level become discrete 
distributions (i.e., histograms). The smearing matrix is 
stochastic: all elements are non-negative and the sum of 
elements in each column is equal to 1. Thus, the matrix 
columns are the probability density functions that relate 
each bin in a histogram at the particle level to the bins 
in a histogram at the reconstruction level. We split the 
full AS {A(j>) range [0, tt] into eight bins and fix two bins 
at small AS* (A0) angles as — 0.7 and 0.7 — 1.2 radians. 
These two bins arc the most sensitive to a contribution 
from DP scattering, and their widths are chosen as a com- 
promise between sensitivity to DP events and the size of 
the relative statistical uncertainty. The sizes of the other 
six bins are varied to minimize the ratio of maximum and 
minimum eigenvalues of the smearing matrix, defined as 
the condition number. This ratio, greater than unity, rep- 
resents the scaling factor for the statistical uncertainties 
arising from the transformation of the differential cross 
sections from the reconstructed to the particle level. The 
same binning is used for the reconstructed and particle 
level distributions. We build the smearing matrix using 
reconstructed and particle level events in the reweighted 
MC sample, described in Sec. IV Bl To decrease the sta- 
tistical uncertainties (^sta.t, ) , we use a Tikhonov regular- 
ization procedure [47H50[ for the matrix. This unfolding 



U 



(O.562) +(<5stat.2) . (6) 



These two bins, being the most sensitive to contri- 
butions from DP scatterings, arc the most important 
for our analysis. We perform the regularization of the 
smearing matrix by adding a non-negative parameter a 
to all diagonal elements of the smearing matrix. The 
matrix columns are then re-normalized to make the ma- 
trix stochastic again (a = is equivalent no regular- 
ization, while for a — >■ 00 the smearing matrix becomes 
the identity matrix). The smallest uncertainties U are 
usually achieved with a = 0.3 — 0.5. An estimate of 
the unfolded distribution is obtained by multiplying the 
histogram (vector) of the measured AS and Acj) distri- 
butions by the inverted regularized smearing matrix. We 
use the sample of the reweighted MC events to get an 
estimate of the statistical uncertainties and the bias in 
bins of the unfolded distribution. To accomplish this, we 
choose a MC subsample with the number of events equal 
to that of the selected data sample and having a discrete 
distribution (histogram) at the reconstruction level that 
is almost identical to that in data. We randomize the MC 
histogram at the reconstruction level repeatedly (100,000 
times) according to a multinomial distribution and mul- 
tiply this histogram by the inverted regularized smearing 
matrix for each perturbation. We obtain a set of un- 
folded distributions at the particle level. Using this set 
and the true distribution for this MC sample, we estimate 
the statistical uncertainty Sstat and the bias b as the RMS 
and the mean of the distribution "(true— unfolded)/true" 
for each AS* (Acfi) bin. The unfolded distribution is then 
corrected for the bias in each bin. We assign half of 
the bias as a systematic uncertainty on this correction. 
The overall unfolding corrections vary up to 60%, being 
largest at the small angles. The total uncertainties, esti- 
mated in each bin i for the AS* and A0 distributions as 

(0.560^ + (4tata)^ vary between 10% and 18%. 
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